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One of the most exciting discoveries in astrophysics in the past few years has been the 
presence of an anisotropy in the cosmic microwave background (CMB) [|T|. It is generally 
believed that this anisotropy is induced by metric fluctuations generated in the early uni- 
verse. CMB anisotropy measurements could thus provide an invaluable test of their origin. 
However, before one can really compare the theoretical prediction of the induced anisotropy 
in any specific model to the experimental data, one has to know within the model the 
theoretical error of the calculation. This error, which is often called cosmic variance, is 
associated with the underlying statistics of the fluctuations. While a theoretical calculation 
gives a definite ensemble-averaged mean value for the anisotropy, a terrestrial measurement 
only measures the anisotropy of a single sample universe. If one uses the experimental data 
to constrain a model, one will then have to take into account the theoretical error of this 
mean value. 

In CMB measurements, the measured temperature anisotropy 6T/T{e) (where e is a unit 
vector pointing to the celestial sphere) is usually expanded in terms of spherical harmon- 
ics, 6T/T = J2i,m(^imYim, from which a two-point temperature correlation function can be 
constructed. After averaging over the celestial sphere (sky-averaging), the function is given 
by 

m ^ (^(ei)^(e2)).., = l^Y.CiPi{cose) , (1) 

where 6 is the separation angle, and Ci = J2m=-i Am'^^rn is the anisotropy power spectrum. 

Estimates of this cosmic variance in the anisotropy power spectrum have been made 
previously They claim that all independent gaussian random variables 

satisfying 

/ t \ {C'l)enSj- (- 

{O'lm^l'm'lens — 2/ -|_ X ' V^J 

where {)ens denotes taking ensemble average. As there are 2/ -|- 1 of the a/^'s for each I, 
using chi-squared statistics one then gets 

{AQr ^ {{Q - {QU^r) = ^i^. (3) 

The variance ACi thereby decreases with l/\/2l + 1. 
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The purpose of this paper is to present a first principle calculation of the theoretical error 
in Ci for the tensorial component of the CMB anisotropy caused by a gravitational wave 
(GW) background. We concentrate on the tensor contribution because after the Planck 
era the GW couples extremely weakly to matter and is essentially a free field propagating 
on a curved metric. During inflation, quantum fluctuations in the gravitational field are 
stretched out of the horizon where they freeze and remain constant in amplitude. When 
they re-enter the horizon much later they appear as classical GW's. These GW's are often 
approximated as a classical stochastic background in the literature [0,^ and Eq. is 
based on this approximation. Using the Heisenberg representation approach [§, we do not 
make this approximation and are able to calculate the cosmic variance explicitly. As ours 
will be a quantum mechanical calculation, we shall use {O) to denote the expectation value 
of an operator. In classical language, this is equivalent to taking the ensemble average. 

We assume a flat Robertson- Walker metric. The quantized graviton field operator hij 
may be written as 

q s=l 

where A is an overall constant whose value is not important for our purposes, r] and R are the 
conformal time and cosmic scale factor respectively while pfj{q) is the circular polarization 
tensor for a GW in state s. a'^^ijf) and a^^irf) are lowering and raising operators in the 
Heisenberg representation with their evolution governed by the hamiltonian P] 

^ = E E + <l<^-i{vy.q{v) + MV) [a^ivy-iiri) - a^{ri)a^_^{ri)\ } , (5) 

q s=l 

where (j{ri) = iR' /2R and the prime denotes derivative with respect to rj. Note that it is 
in a simple quadratic Gaussian form. The Heisenberg evolution equations, = [a|, if], 
ia'^y = —[ap,H], can be solved exactly using the Bogolubov transformation: 

afiv) = u'q{v)a^Hvo) + vl{v)a%ivo) , (6) 

where rjo is some initial time. Then m^' = —iqu^ + v^R'/R, = —iqVg + u^R' /R, with the 
initial conditions Ug(?7o) = ^g('7o) =0 0. 



In the Schroedinger approach found in |]T^ the wave-function itself is time dependent 
while the operators are time independent with the field operator written as 



R 



(7) 



To find the relationship between h and the u-v coefficients, we use Eq. in Eq. (|^) and 
compare with Eq. (|^), giving = Ug + Vg. Hence, h^" + (g^ — R" / R)h^g = and the initial 
conditions for and give h'^qirjo) = 1 and h^g'^rjo) = —iq. 

Following 0, the classical temperature fluctuation 6T/T of the CMB is now replaced by 
a field operator, 6T/T = — /J" dKh[j{rj, x)e*e-'/2, where the lower (upper) limit of integration 
represents the point of emission (reception) of the photon. Note that aim, Ci, and C{9) 
defined in Eq. (|l]) now become field operators with 
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{C{e)) = (0|^(e0^(e2)|0) = - / rfAi / dA, e\^efel^ei^ 



d d 



■ (0 1 (r^i , f 1 ) (r?2 , ^2 ) 1 0) 



dr]i drj2 

where |0) is the vacuum state of the system at time tiq. Using Eq. 

.(2) 



(8) 
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«1J1,*2J2(^1'^2) = {^\Khi.^l)Kj2i.^2)\Q) 



Rivi) 



we find 



<(^2) +v'Jr]2 



Rim) 



(9) 



and {C{6)) is dependent on the two point Green's function for the field. 



Notice that G) 

.(2) 



(2) 



3^1,3^2) = ^j2i2,nii('''2 5 a;i). For a classical field, however, one would 



expect Gj^iji,i2i2(^i'^2) = G\l]^^i^j^{x2,xi). From Eq. (g), this holds only when + is 
real. Of course, h^g = u^g + Vg will not, in general, be real given the initial conditions for 
Ug + Vg. It can be shown, however, that if immediately after inflation there is a transition 
to the radiation-dominated era, hg is essentially real for all relevent GW wavelengths up 



to an overall g^dependent constant phase factor after inflation |jTT],|T2[. This phase can 
be considered as a random spatial phase which determines the location of the node of 
a Fourier mode of the wave. It is randomly distributed since quantum fluctuations during 



inflation assign equal probability to modes which differ only by a spatial translation |]13|, P 



Consequently, the usual "scale-invariant" initial condition for a GW generated from inflation 



at the end of the inflationary era at rj = 0, (up to an irrelevant overall constant phase) is 
given hj h^g/R = 1, (h^g/R)' = as — > 0. Subsequent temporal evolution of the wave can 
be well approximated by solving the equation of motion for /i^ under these initial conditions 
p!5| . More importantly, the use of classical GW's after inflation in the calculation of (Ci) is 



justifled, and (Ci) is determined only by the temporal phase of the GW. 
Using the identity from p, 

J dmU^)e''-%{q}e'e^ = Mq\x\)H'm , (10) 
where H^, = (5+ - ^5-)^ + i^t + ^^X, and 



1/2 



25r _^ 5r+2 
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we obtain 



(2/ - 1)(2/ + 3)(/ + 1/2) (2/ + 1)(2/ + 3)(/ + 5/2) 



{Ci) = IZ(aL«im) = / dAi dA2^^gi{vi,V2) , (12) 



(2/- 1)(2/+ l)(/-3/2) 

(11) 



where 



9i(Vi,V2)=A^EMQ\mMQ\x2\)^^. (13) 

Note that in Eq. we have taken h^^ as a real function obeying the scale invariant initial 
condition. 

Next, to calculate the variance in Ci we note that (Cf) depends on the four-point Green's 
function for the fleld. Using Eq. (|), 

^iiji,...,j4j4 (-^l' • • • 5 -^4) = ^iiii,j2i2 ("^l' '^2)Gjjjj^ j^j^ (X3, X4) + 

^ilil,j3i3(^l' ^3)^i2j2Mj4(^'^^ ^4) + 

^ilil,j4i4(^l'^4)G-2]2^i3j3(X2,X3) . (14) 

This form of factorization is expected from a Gaussian theory. From this, it is then easy to 
show that 

{iCi-{Q)r) = 2j:\{aLai^')\'. (15) 



From Eq. ([T0|), it is easy to show that J^sH^H^, = 2Smm', indicating that each component 
in m is orthogonal to each other. In addition, due to the isotropy of the GW background, 
equal {al^aim)'s for a fixed / would be expected. Hence, it follows from Eqs. (|12D and (|15D 
that 

(AC,)^ - (16) 

which is exactly the classical result (H). 

In conclusion, we have presented a calculation of the theoretical error in Ci based on 
quantum field theory. We found from the isotropy and Gaussianity of the theory that the 
fractional error in Ci decreases with 1 / \/2l + 1, as expected in the classical approximation. 
We have, however, considered only the tensor contribution to the anisotropy of the CMB 
in this paper. As it is believed that the main contribution to the CMB anisotropy is most 
probably due to the scalar perturbations, a calculation of the cosmic variance for the scalar 
part is needed before any definitive characterization of the experimental data can be done 
with any certainty. 
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